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We analyze tunneling of Cooper pairs across voltage biased asymmetric single-Cooper-pair transis- 
tors. Also tunneling of Cooper pairs across two capacitively coupled Cooper-pair boxes is considered, 
when the capacitive coupling and Cooper pair tunneling are provided by a small Josephson junc- 
tion between the islands. The theoretical analysis is done at subgap voltages, where the current- 
voltage characteristics depend strongly on the macroscopic eigenstates of the island(s) and their 
coupling to the dissipative environment. As the environment we use an impedance which satisfies 
Re[^((j;)] <C Rq and a few LC-oscillators in series with Z{ijj). The numerically calculated I — V 
curves are compared with experiments where the quantum states of mesoscopic SQUIDs are probed 
with inelastic Cooper pair tunneling. The main features of the observed I — V data are reproduced. 
Especially, we find traces of band structure in the higher excited states of the Cooper-pair boxes 
as well as traces of multiphoton processes between two Cooper-pair boxes in the regime of large 
Josephson coupling E,j ^ Ec. 

PACS numbers: 74.50.-|-r, 73.23.Hk 



I. INTRODUCTION 

A voltage biased small Josephson junction (JJ) has 
been shown to be a good probe of mesoscopic physics. In 
recent years it has been used, for example, in the detec- 
tion of resonances in the electromagnetic environment^'^ 
and noise spectroscopji^. The theory of inelastic tun- 
neling, known as the P(£^)-theory, describes I — V char- 
acteristics resulting from incoherent tunneling of Cooper 
pairs, or quasiparticles, across the small JJ and simul- 
taneous energy exchange between the tunneling particle 
and its electromagnetic environment, which is described 
by a set of LC-oscillators. The standard P(£')-theory 
cannot, however, be used in the case of a non-Gaussian or 
anharmonic environment. In this paper, a suitable model 
will be constructed to account for the anharmonicity of 
the environment. 

This paper gives a quantum description for a system 
which is designed to probe the excited states of a Cooper- 
pair box (CPB), or coupled boxes, by a small J J. We 
model the quantum evolution of a voltage biased asym- 
metric single-Cooper-pair transistor (SCPT) or a circuit 
consisting of three JJs in series with a small middle J J. 
The idea is, as in the P(i5)-theory, that the small JJ 
is probing the eigenstates of the CPBs, which are then 
seen as current peaks at certain voltages. This is possible 
since under the voltage bias well above the supercurrent 
peak, but still at subgap region, the tunneling of a single 
Cooper pair across the small JJ is possible (nonvirtually) 
only if the environment is able to absorb the energy 2eV 
released in the tunneling. 

The environment of the small JJ consist of a CPB and 
a continuous spectrum of LC-oscillators describing dis- 
sipative quantum mechanics induced by high frequency 
resistive properties of the leads and possible spurious res- 
onators in the transmission line or materials nearby the 



island. In resonant situations the dynamics involve both 
excitation and relaxation of the CPB eigenstates and one 
is, in principle, able to get information of both the ener- 
gies as well as the relaxation times of the excited states. 

Experimentally, the spectroscopy of the eigenstates us- 
ing a small J J as a probe have been done by Lindell et. 
al. and the results are reported in Refs. IsIgI Fti In this 
set of experiments, traces of excited states, their anhar- 
monicity and expected band structure were found from 
the measured I — V characteristics. However, several 
unexplained phenomena seen there were the main moti- 
vations for writing this more detailed description for the 
system. We show that indeed the main features of the 
I — V data can be explained by the quantum mechanics 
of asymmetric SCPTs or coupled CPBs. The model ex- 
plains, for example, the widening of the I ^V resonances 
as result of a band structure of (coupled) CPBs and non- 
constant peak splitting in the experiment of Ref. |a as a 
result of multiphoton transitions between eigenstates of 
two CPBs. 

The paper is organized as follows. In section ^ we 
build a theory describing inelastic tunneling across the 
small JJ when it has an anharmonic element, i. e. a CPB, 
in its environment. In section llTTI we discuss effects caused 
by slow relaxation and section IIVI is devoted to a quanti- 
tative discussion of the I —V characteristics in the case 
of a two CPB environment. Comparison between numer- 
ical calculations and experiments is presented in section 
IVland conclusions arc given in section IVll 



II. 



INCOHERENT TUNNELING OF COOPER 
PAIRS ACROSS ASYMMETRIC SCPT 



We model an asymmetric single-Cooper-pair transis- 
tor by taking the Josephson coupling across the probe 
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this limit the effect of the operator i^EE + -f^int for the 
CPB can be described by a transformation V ^ V + Vf 
in Eq. (O, where Vf describes fluctuations from the av- 
erage value ViiLiS. Therefore one can use an effective 
Hamiltonian for the CPB 
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— HqPB — QintVf, 
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FIG. 1: A voltage biased SCPT. We study the case Eji > 
Ej2 and the smaller JJ is called the probe. The supercon- 
ducting leads (and the spurious environment of the island) 
are modelled through an impedance Z{ij). 



junction into account perturbatively. The treatment de- 
scribes incoherent tunneling of Cooper pairs across the 
small JJ and simultaneous energy exchange between the 
tunneling Cooper pair, CPB and the dissipative environ- 
ment. It is valid if the relaxation rates of the excited 
states are higher than the excitation rates induced by 
incoherent tunneling. 

A voltage biased SCPT is shown in Fig. ^ We are in- 
terested in the case of strong asymmetry, Eji ^ Ej2. 
The charging energy of the island, defined as Eq — 



where Qint — C2Q/CS (we have assumed that Co ^ Cs). 
The autocorrelation function of the fluctuating voltage 
is related to the dissipative properties of the impedance 
Z{ll!) via the quantum fluctuation-dissipation theorem 
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(4) 



Since fluctuations are only a small perturbation to the 
CPB, their effect is to induce transitions between the 
unperturbed states, i. e. the eigenstates of the CPB. The 
transition rate between the eigenstates \i) and |/) is ob- 
tained by the golden rule calculation 
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e /2Ce where C^ ^ Co + Ci + C2, can however have ^^^^^ ^ = E, - Ef ^ fuv is the difference between the 



an arbitrary value. Disregarding the Cooper pair and 
quasiparticle tunneling across the probe (the former is 
taken into account as a perturbation), the Hamiltonian 
of the environment of the probe can be written formally 
as 



if 



corresponding eigenenergies. 

We proceed by noting that the rates (0 define the 
lifetimes, and therefore also the linewidths, of the energy 
levels in the Cooper-pair box. The full width at half 
maximum (FWHM) Aq of the state \a) is then 



Hp, 



HqPB + ^EE + Hint , 



where the CPB Hamiltonian i; 

\2 
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Hqpb — 



(g + Q o) 

2Cy 



Eji cos{ipi), 



(1) 



(2) 



Aa = ft^7/<_a. 



(6) 
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and Qo ~ CoU — (C2 + Co/2)V is the quasicharge, Q the 
island charge, a conjugated variable to the phase differ- 
ence ipi. The Hamiltonian TJee rnodels dissipative elec- 
tromagnetic environment and possible spurious fluctua- 
tors in the system. Its characteristics are fully described 
by an impedance Z{uj), so formally it consist of an infinite 
number of harmonic oscillators. The interaction term 
Hint describes linear coupling between the CPB and the 
electromagnetic environment leading to dissipative quan- 
tum mechanics^*^. 

To account for the anharmonicity and the band struc- 
ture of the CPB, we will proceed slightly differently than 
is done in the P(£')-theory. The idea is to take the anhar- 
monic parts into account separately. This is possible in 
the limit Re[Z(w)] ^ Rq, where the interaction between 
the CPB and the dissipative environment is weak and one 
can use the Born-Markov approximation when describing 
the evolution of the CPB under the Hamiltonian (Q . In 



and the density of the excited states broadens from the 
sum of delta- functions to sum of Lorentzians, i. e. the 
density of states changes as 



E^(^-^")-E^ 
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TT 4(^„ - S)2 + A2 



(7) 



Finally, we include the Josephson coupling Ej2 cos((/32) 
describing tunneling of Cooper pairs across the probe 
junction and simultaneous excitations of its environment 
as (another) perturbation. Using (fi^ — (pi + ip2 + Pz = 
2eVt/h, where (pz is the phase difference across the 
impedance, one finds the time dependent perturbations 
for the positive and negative direction tunneling M± = 
Ej2 exp [±i {ifii +ipz- 2eVt/h)]/2. 

The transition rates due to perturbations M± between 
the eigenstates of the iJenv are effectively described by 
transition rates between the states |a)|env), which con- 
sist of two independent parts: broadened CPB states \a) 
and continuous distribution of the environmental LC- 
oscillators, whose free evolution is described by -ffEE 



(Born approximation). The LC-environment can be 
traced out similarly as in the P(i?)-theory and the golden 
rule transition rates between the CPB states \i) and |/) 
become then 
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A, + A 
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A{Es~Ea~E'Y + {^, + ^fY' 

(8) 

where the P(i?)-function is the same as for a system con- 
sisting of a probe junction with a capacitance C12 = 
(1/Ci + 1/C2)~^ in series with the impedance Z{u)). 
If Z{ijj) is a constant R (^ Rq), the main contribu- 
tion of the P(i?)-fmiction becomes from low energies 
where it is approximately a Lorentzian with a linewidth 
Aonv = 4:TTkBTR/RQ centered at £; = 0^. Therefore it 
convolutes the original transition rates (O to 



r 



A total 



(9) 



i{Ef -E,T 2eF)2 + {A^f^^y ' 



where A'^*''' = Ai + Af + Aonv We see that there are 
two sources of broadening of the resonances: widening 
due to finite lifetimes of the CPB-eigenstates (AqIs) and 
widening due to low frequency fluctuations of the LC- 
environment (Aonv)- The separation of the CPB and its 
environment holds also for the case where the environ- 
ment has several modes which are nondegenerate with 
the CPB eigenstates. The degeneration, or almost de- 
generation, would lead to similar splitting of the states 
as described in section Hvl 

Each transition will occur simultaneously with a trans- 
fer of 2e of charge across the system and the current is 
therefore 



2e 






r 



/*- 



(10) 



where the probabilities Pi for occupancies of the CPB 
eigenstates are given by the canonical equilibrium distri- 
bution. If ksT <t:Ei- Eo and Z{lu) = R then 



/(y)=2e^r;^„ = 5:^|(/|e- 
/ / 

A total 

fV 

4{Ef -Eo- 2eVf + {A^^y ' 



(11) 



One sees that I —V peaks can be identified with energy 
levels of the environmen'feSii^, which in this case is the 
Cooper-pair box. 

We have verified that the I—V characteristics obtained 
from Eq. lfTT|l reduce to the ones obtained from the P{E)- 
theory, if the larger JJ is described as an LC-oscillator. 
However, Eq. \i\\ is also valid for JJ-environment with 
evident anharmonicity or band structure, and therefore 
is not limited to the harmonic approximation. 



III. EFFECTS DUE TO SLOW RELAXATION 

For the expression @ to hold, it is vital that the sys- 
tem relaxes rapidly to the ground state, since the golden 
rule calculation is justified only if the excitation rates of 
the CPB eigenstates, induced by the probe, are smaller 
than the relaxation times, caused by the CPB's coupling 
to the dissipative environment. The irreversible interac- 
tion with the dissipative environment has to "cut" the 
evolution to the excited state quickly, otherwise the tun- 
neling across the probe would turn from incoherent to 
coherent. On the other hand, in the opposite case of 
very slow relaxation, one would obtain Rabi oscillations 
between the CPB eigenstates |0) and |a) when initially 
starting from the state |0) with 2eV = Ea — Eq. This 
limit can also be analyzed in the Born-Markov approx- 
imationi^, but generally the problem needs an analysis 
of the time evolution of the whole density matrix and 
Markovian approximation cannot be usedi^. 

To obtain approximative results in all regions, we use 
the model derived in section^ with modified probabili- 
ties for occupations. We redefine the diagonal elements 
of the density matrix by the ones obtained from the equi- 
librium master equation 






+ 7/^.)] = 0, (12) 



for each i. The 7:s are the relaxation rates caused by 
the fluctuating voltage across the CPB, Eq. (O, whereas 
the r:s are the rates induced by Cooper pair tunneling 
across the probe, Eq. ^. The method assumes that all 
the sequential transitions are independent of each other, 
which is not always true. However, the method reduces 
to the one considered in section ^ when the relaxation 
dominates the excitation and, according to our numerical 
calculations, gives similar results for the first order tun- 
neling processes (single Cooper pair tunnels across the 
probe with simultaneous excitation of the CPB) even 
in the regime of very slow relaxation, as long as the 
SCPT is highly asymmetric. Therefore it is safe to as- 
sume that the first order processes are well approximated 
by Eq. IjlUI) with the equilibrium probabilities obtained 
from Eqs. (|12|l . Also, there is no experimental evidence 
of higher order resonances, which might be due to their 
weak nature to be washed out by the so called Zeno ef- 
fecO^. 



IV. TWO CAPACITIVELY COUPLED CPBS 

To generalize the treatment of the preceding sections, 
we do a perturbative treatment for three JJs in series, 
where the middle one acts as a probe. The configura- 
tion can be seen to consist of two capacitively coupled 
Cooper pair boxesi^, where the capacitive coupling is in 
parallel with a small tunneling element. Since we use 
similar models for analyzing the experiments in section 
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FIG. 2: Three Josephson junctions in series with the voltage 
source. One of the large JJs could as well as be an LC- 
oscillator, modelling the effect of a spurious resonance. 



IVI we concentrate on the characteristics of this model a 
bit deeper. We also note that one of the larger JJs could 
as well be an LC-oscillator describing spurious resonance 
at frequency ujp — l/vLC in the environment. If the en- 
ergy quantum tnop is almost the same as any excitation 
energy En ~ E,n between two relevant eigenstates of the 
CPB, or the state is long living, it cannot be modelled by 
the P(£^)-function in Eq. 0, but the following treatment 
is valid. 

The system in consideration consists of three JJs in 
series connection with the voltage source and the smallest 
junction is in the middle, Fig. 13 Following the steps 
done in section^ we first neglect the Josephson coupling 
energy of the probe and write down the Hamiltonian of 
two capacitively coupled Cooper pair boxes 



ff2CPB — 



{Ql + Qo) , (Qr + Qo) 



2Cl 

Eji cos{ipi) ~ Ej3 cos(<^3; 



2Cr 

Ci23QlQr 

C1C3 ' 



(13) 



where ipi and Ql are conjugated variables and simi- 
larly for (/33 and Q^. The capacitances of the islands 



where Cgi is the ca- 



are (assuming that Cgi <^ C, 
pacitance of the gate i) Cl = Ci + C23 where C23 — 
(I/C2 + l/C3)-\ Cr ^C3 + C12 and C123 = (1/Ci + 
I/C2 -I- l/Cs)"^. The quasicharges become then Qq w 
CgiUgi + Cg2Ug2C2/CR " ^2 F and Q'o « -Cg2Ug2 - 
CgiUgiC2/CL — C2V. One sees that the Hamiltonian 
consist of two CPB Hamiltonians, which are then coupled 
by the last term in Eq. Ijl^l) . Similarly as in section ITU 
we first calculate the linewidths of the eigenstates of the 
coupled system (but now Qi^t = QlC2/Cl + QrC2/Cr), 
then use the fact that the phase difference ipj: is a classical 
variable and take the tunneling across the probe into ac- 
count perturbatively by considering the broadened states 
of the coupled CPBs and the environmental states sepa- 
rately. The current is obtained from Eq. UlUfl. similarly 
as before. 

The behaviour of the eigenstates and energies of the 
Hamiltonian H13|l can be analyzed analytically in the 
limit Eji,Ej3 > e'^/2CL,e'^/2CR. For simplicity let us 
assume that Eji = Ej3 = Ej and Cl = Cr « Ci = 



C3 = C. Two "splitting" effects contribute to the fi- 
nal energy level structure of the coupled CPBs. First, 
in the harmonic approximation of Eq. (|13|) the CPB:s 
behave as LC-oscillators. The degeneracy of the iden- 
tical LC-oscillators is removed by the interaction term 
Ci23QlQr/CiC3 w C2QlQr/C^ (presuming that C2 < 
C) , and diagonalising the quadratic Hamiltonian one sees 
that the system behaves as it would consist of two inde- 
pendent oscillators with the original inductances but with 
capacitances C± — C'^/{C ± C2). For small C2/C this 
leads to mode frequencies u>pzLujpC2/2C. Secondly, if one 
takes into account the first nonharmonic terms —EjLpf/Al 
in the Hamiltonian H13|l and sets C2QlQr/C^ — 0, one 
obtains also energy level splitting effects due to com- 
bined energy levels of two anharmonic oscillators. The 
energy levels of single CPBs become Ei = fwjp — Eq, 



E2 



2hujp — 3Ec, E3 — Shwp — 6Ec ■ ■ 



New levels 



appear due to simultaneous excited states of the boxes 



12*) = 11)11) 



|3*) = -^(|1)|2) + |2)|1)) 



(14) 
(15) 



with the corresponding eigenenergies £'2* — 2hujp — 2Ec 
and E3* = Shuip — AEc- The energy level 2?iajp therefore 
"splits" into two nearby energy levels £'2 and £2* ■ 

When both of the above effects are included, more mix- 
ing of the states is obtained. The n:th excited state splits 
into 71-1-1 states and, for example, the state |1) splits into 
states (using the first order perturbation theory) 



|l.) = i=(|l)|0) + |0>|l» 
|la) = i=(|l)|0)-|0)|l)) 



(16) 
(17) 



with the eigenenergies Eis = hujp — z — Ec and Eia 



hujp + z - 
state 1 2) 



Eq, where z = C2h^pl2C. Similarly for the 



2.)=ci(|2)|0) + |0)|2)+c+|l)|l)) 
^(|2)|0)-|0)|2)) 



\2a 
|2*)=C2(|2)|0) 



|0)|2)+c-|l)|l)). 



(18) 
(19) 
(20) 



where i?2s 
E2' = 2hw. 



2huJr, 



5Ec/2 



E2 



2hu!„ 



■5Ec/2 + z', 



VW 



Z', ii2a = '^n^p - 3Ec, 

i-Ec/2±z')/V2z,z' = 

fj -r 16z^/2 and Ci are normalizing factors. These 
states give both behaviours discussed in the preceding 
paragraph as the limiting cases of Ec -^ and z ^ 0, 
respectively. The states |la) and \2a) do not lead to cur- 
rent peaks since they contain antisymmetric combination 
of the states and therefore the elements | (/| exp{±i{ipL + 
(fR,))\0)\ vanish. We still pick up the energies of the 
states £"3.; ~ Shujr, 
£3* 



5£c - z - 



^/W 



4z2 _ 2Ecz and 



Shujp - 5£c - z+ y/E^ + 4^2 - 2Ecz. 
For higher anharmonicity (£j ~ Ec), where the band 
structure will become evident for the eigenstates and 



the perturbative treatment of the cosine-potential is no 
longer valid, we have to resort to numerical solution. We 
calculate the eigenstates by diagonalising the Hamilto- 
nian in a product basis of two noninteracting CPBs, for 
given values of gate voltages (quasicharges). The low- 
est eigenstates can be obtained quite accurately from an 
economical sized matrix equation, since the eigenstates 
are usually close to the states of this basis (due to small 
capacitive coupling), justifying also the "product state 
labelling" of the final states. After the Hamiltonian is 
diagonalised (and the transition rates have been calcu- 
lated) one has to solve Eqs. H1U|I and (|12|l for each value 
of V to obtain the I — V characteristics for given Ugi . 

In section IV Bl we model experiments using a similar 
circuit but including also two extra iC-oscillators in se- 
ries with the three J J system. The previous I — V charac- 
teristics of the coupled CPBs are still preserved but mul- 
tiphoton transitions with the external LC-oscillators are 
also obtained, which is the motivation for this procedure. 
In practice, the extra LC-oscillators can be modelled as 
J Js and the Hamiltonian of the system can be written as 



H = ^^{C ^)kiqkqi-^E,ncos{Lpi) 



(21) 



k,l 



where q^ = Qk — Qk+i, Qk is the charge gone through 
the fc:th JJ, a conjugate variable to ipi, C is a capacitance 
matri^i. This Hamiltonian fully determines the energy 
bands, i. e. the ranges where peaks can occur in / — 
V characteristics. In order to determine further details 
such as the peak positions for given gate voltages, one 
has to complete the Hamiltonian with linear terms in 
charge. Such terms result from voltage sources and can 
be deduced from single tunneling eventsi^i^. However, 
it turns out that in the experiments to be analyzed, the 
quasicharge is averaged over all values and the resonance 
positions become immune to these terms. 

The relaxation rates due to photon emission to 
the electromagnetic environment Z{lo) are also deter- 
mined by linear terms, through the fluctuating operator 
QintVf = J2i'^iQi^f- For the LC-oscillators this relax- 
ation channel does not lead to observed rates and must 
be enhanced by introducing resistors in parallel with L 
and C. An analogous procedure is to take the coeffi- 
cients qlc as fitting parameters and use the operator 
J2i O'LaQLCiVLCi S.S & perturbation, where the fluctua- 
tions Vlci are uncorrelated but have the original prop- 
erties. For the large JJs a^'s are theoretically defined 
by the ratios C2/Ci/jj, as was seen in the beginning of 
this section. The ratios are locked when fitting the ob- 
served energy level structure of the coupled system. But 
in real systems the capacitive coupling of CPBs can be 
effectively reduced by decoherence effects, for example 
by thermal fluctuations or dissipation, leading to a de- 
crease in the "observed" C2 . Also effects related to mate- 
rials nearby the CPB islands seem to be able to increase 
relaxatioR^^. Therefore in modelling the relaxation rates 
of coupled CPBs, one is forced to take the corresponding 
coefficients a^ as independent fitting parameters. 



V. COMPARISON TO EXPERIMENTS 

The I — V characteristics of similar systems as dis- 
cussed above were measured experimentally by Lindell et. 
al. and reported in Refs.|^O0 In these experiments, dif- 
ferent kind of environments for the probe junction, con- 
sisting of one or several SQUIDs and two or four leads, 
were used under different magnetic fields and gate volt- 
ages. Using SQUIDs as the large JJs, the system could be 
studied in situ from the harmonic behaviour {Eji ^ Ec) 
to the region where the anharmonicity and band struc- 
ture become crucial {Eji ~ Ec), by applying magnetic 
flux to the SQUID loops. This property also helped in 
the analysis of the data, since resonances coming from the 
spurious environment did not react to the applied mag- 
netic field, at least not in the same way as the resonances 
coming from the CPBs. 

The probe junction current as a function of voltage and 
external flux is shown in Fig.|21as a 2D-surface plot. The 
dominant current peaks show periodicity as a function of 
the flux $ through the SQUID loop with the period of the 
flux quantum <I>o = h/2e. This clearly points that they 
are originating from the SQUIDs and allows the identifi- 
cation of the different SQUID excitation states from the 
more complicated I — V structure due to the rest of the 
electromagnetic environment. In addition to the peri- 
odic structures due to the SQUID environment, one can 
see additional states with longer and non-constant peri- 
ods. It is likely that these are due to large, additional, 
Josephson junctions that are created in the two-angle 
evaporation technique used to fabricate the sample. The 
patterns have the Fraunhofer/Airy behavior as expected 
for a large Josephson junction that is penetrated by a 
magnetic fieldi^. 



A. The 1-SQUID Experiment 

We begin by studying the sample that has the con- 
figuration of the asymmetric SCPT (Fig. ^. The first 
thing in the fitting procedure is to identify which of the 
resonances in the I — V characteristics are coming from 
the CPB eigenstates, which from spurious fluctuators 
and which from simultaneous excitations of both. The 
I — V peaks in this " 1-SQUID experiment" consist of a 
set of flux dependent double peaks and several static res- 
onances, from which the most important is at Vlc ~ 13 
^V, i. e. ljo/2tt w 6.3 GHz, see Fig.0] Its second excited 
state is seen at 2Vlc ~ 26 fiY (not shown in Fig. ^J 
and therefore it is not a 2-state fluctuator. The reso- 
nance is important since it explains the double structure 
of the first two flux-dependent double peaks: the lower 
resonance of each double peak is due to tunneling of a 
Cooper pair across the probe and simultaneous excitation 
of the CPB and the higher resonance of each double peak 
is due to tunneling and simultaneous excitations of the 
CPB and the LC- resonance (a multiphoton transition). 
Two observations support this idea. First, the peak split- 
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FIG. 3: A 2D-surface plot of the measured probe junction 
current as a function of voltage and flux through the SQUID 
loop. The resonances originating from the SQUIDs show 
$o-periodicity. The figure also shows resonances of Fraun- 
hofer/Airy type with weaker dependence on <I>, and multi- 
photon transitions. The I — V characteristics of this 4-SQUID 
sample are analyzed more quantitatively in section IV Bl 



ting is constant as a function oi Eji{^) — Eji cos($/$o) 
and this constant is the same for the first and the second 
double peaks and also equals Vlc, reflecting the same 
excitation energy difference 2eVLc, see Fig. 0] Secondly, 
when compared to the first peak of a double peak, the 
relative area of the second one is approximately the same 
for the first two double peaks, indicating a similar extra 
factor (matrix element) related to the latter resonance of 
a double peak. The theoretical model used in fitting is 
the same as in section IIVI except that one of the J Js is 
replaced by an LC-oscillator. 

The third double peak is not, however, consistent with 
similar interpretation, since the peak splitting is smaller 
than the previous ones at <& = and increases with in- 
creasing $. Also, the two resonances of the double peak 
have almost equal areas, but the multiphoton transition 
to the external LC-oscillator should have almost vanish- 
ing area and should not even be seen. The peak splitting 
can be explained by the band structure of the CPB's 
third excited state, assuming that the observed I — V 
curves are certain averages of the quasicharge-space and 
band edges are highlighted due to van-Hove-like singu- 
larities. The observed splitting indeed follows the reso- 
nances obtained from band edges, as seen in Fig.^] Still, 
the physical reason for the "escape" of the quasicharge 
is unknown. No gate dependence for the positions of the 
I — V peaks is seen when Eji > Ec, supporting the idea 
of the "running" polarization charge. It is, interestingly, 
returned in the limit Eji < Ec at higher voltages as 
charging effects, when the gate-dependence of the ground 
state energy becomes observable. 

Indications of the band structure of |2,0) and even 
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FIG. 4: The positions of the resonances in the 1-SQUID 
experiment (data points) compared with those resulting from 
the model of a SQUID and an LC-oscillator in series with 
the probe junction (lines). The resonances have been labelled 
by the corresponding transitions behind them. For example, 
|1, 0) means that the resonance occurs due the transition from 
the ground state to the product state, in which the CPB is 
in its first excited state and the LC-oscillator stays in its 
ground state. The energy levels of the CPBs are actually 
bands (shaded), and the resonances due to the band edges 
q = 0,q — e have been plotted explicitly. The parameters in 
the numerical modelling are summarized in table Q] 



1 1, 0) excitations are seen in the experiment for <i> > 0.3<i>o 
and <i> > 0.4$oj respectively, but instead of clear split- 
ting the resonances widen and become fluctuating. This 
broadening is also expected theoretically, as shown in 
Fig. 21 and looks like to be caused by random quasicharge 
fluctuations. Unfortunately the noise induced by the en- 
vironment exceeds these current peaks for $ > 0.35<i>o 
and therefore no clear evidence of these bands is ob- 
tained. 



Finally, we note that the linewidths of the resonances 
are similar for the lowest resonances ~ 8 fj,Y, which can 
be fitted using the values Rq = 100 ft and T = 0.4 K 
(leading to Acnv/2e ~ 3.5 /iV) and the independently 
measured value Ej2 = 8.5 /leV. The effective tempera- 
ture is quite high since the temperature of the environ- 
ment was ~ 0.1 K. The model indicates that the lowest 
resonances are in the slow relaxation regime, which is 
consistent with the observation that the maximum cur- 
rent of the 1 1 , 0) resonance decreases when the magnetic 
fiux is increased, see Fig. |S1 
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FIG. 5: The maximum current of the resonance |1, 0) in the 
1-SQUID experiment (triangles) compared with that resulting 
from the model of a SQUID and an LC-oscillator in series with 
the probe junction. Different from Fig. 2] the LC-oscillator 
( "I/C2" in table describes a small I — V peak seen at a 
voltage « 90 \iSf (not shown in Fig. UJ. In the region "I> > 
0.3$o the state |1,0) suffers from a slow relaxation and the 
current is determined by the photon emission to the ohmic low 
frequency environment. A resonance occurs nearby ~ 0.2$o 
when the states |1, 0) and |0, 1) are coupled, the latter having 
a faster relaxation to its ground state. The current between 
$ = and <1> = 0.1 $0 seems also to be enhanced, probably 
due to further entanglement with the environment. In the 
calculation of the relaxation rates we have used a perturbation 
(Ca/COOsQumT// + 0.2Q£cV7' (sectionEVJ. 
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FIG. 6: Left: Schematic drawing of the 4-SQUID experiment. 
The SQUIDs are drawn as JJs. In this situation, the two 
SQUIDs (at the same side) behave as a single JJ, but with 
double the coupling energy Ej and capacitance C compared 
to the individual SQUIDs. Right: The resulting theoretical 
model of the system. 



B. The 4-SQUID Experiment 

The second sample to be studied consists of four leads, 
four SQUIDs and a probe junction, see left side of Fig.El 
Since the two SQUIDs on the same side of the probe 



behave as a SCPT and the phase difference 9 across this 
component relaxes to the minimum energy value 6 — Q 
(its classical dynamics is highly damped and no bias is 
present) , the two SQUIDs behave as a single J J and one 
arrives at an equivalent circuit of two JJs and a probe in 
series connection; the model discussed in section HVI 

Also for this sample, the I—V characteristics consist of 
flux dependent double peaks and a few static resonances. 
The noise of the environment in the limit Ej\ ^ Ec is, 
however, much smaller than in the 1-SQUID sample and 
therefore a more detailed analysis can be done. Again, 
the first static resonance is seen at Vlc ~ H mV i- e. 
cjo/Stt « 5.3 GHz, see Fig. 13 and its second excited state 
is seen at 2Vlc (not shown in Fig. [Tj). We include also a 
resonance seen at Vr « 123 fiV to the model as another 
LC-oscillator. The first flux dependent double peak can 
again be explained a as a plain excitation of the CPB, and 
a multiexcitation of the CPB and the (smaller frequency) 
LC-oscillator. The second double peak, however, is not 
consistent with this assumption since the peak splitting 
is much smaller than UlCj 8.5 /iV, and the relative areas 
of the peaks, 1 : 0.4, differ essentially from the ones ob- 
served for the first double peak, 1 : 0.13. Instead, a bet- 
ter explanation for the second double peak is the energy 
level structure of two coupled CPBs in the anharmonic 
region (section llV|) : the peak sphtting occurs due to reso- 
nances of the |2, 0, 0) and |2*, 0, 0) excitations, and using 
the parameters given in table Q] one obtains a peak split- 
ting 8.5 //V and an area ratio 1 : 0.25. The resonances 
corresponding to the multiphoton transitions |2, 1, 0) and 
1 2*, 1,0) are also seen as weak peaks in agreement with 
the model. The third double peak splitting is then au- 
tomatically explained as excitations to the states |3, 0, 0) 
and |3*, 0, 0), giving the peak splitting 13.2 /zeV and area 
ratio 1 : 1. The experimental values are 15.0 fxeV and 
1:1. Note, that in this sample the band structure of the 
state |3,0, 0) is neglible at $ = 0, and therefore it does 
not explain the third double peak. The multiexcitation 
|1,0, 1) is also seen in Figs. 01 and [7| justifying generally 
the multiphoton interpretation. 

We conclude that the area ratio and the peak splitting 
of the first double peak can be fitted by changing the 
properties of the external LC-oscillator but at the same 
time the resonance at F = Vlc has to be fitted also, 
whereas the other splittings and areas are determined by 
the charging energy Ec of the island(s) and the capaci- 
tance of the probe junction C2 (section Hvj) . Also adding 
a slight asymmetry between the two SQUIDs can "fine 
tune" these values and we have used here a 1% difference. 
From Fig.|Slone can see that, not only the calculated po- 
sitions of the resonance peaks, but also the corresponding 
areas are quite similar for the model and experiment for 
this choice of fitting parameters. 

It is interesting to study the effect of band structure in 
this experiment. Again, a change in the gate voltage did 
not result in a change of the resonance positions, even 
in the region Eji ~ Ec- Instead, the resonances origi- 
nating from SQUIDs widened and changed from smooth 
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FIG. 7: The positions of the resonances in the 4-SQUID 
experiment (data points) compared with the edges of the en- 
ergy bands (lines) and the band structure (shaded) calculated 
from the model of two SQUIDs and two LC-oscillators in se- 
ries with the probe junction. The experimental data is based 
on I — V curves of which two examples are shown in Fig. |H] 
Also shown is a schematic diagram of the model circuit and 
the code used for labelling the states. The parameters in the 
numerical model are summarized in table |T| 



Sample 


i5,7i(MeV) 


Ej2{^ieV) 


Ci(fF) 


C2(fF) r(K) 


1-SQUID 
4-SQUID 


390 (188) 
483 (544) 


8.5 (8.5) 

3.6 (3.6) 


4.6 (5.7) 
6.45 


0.7 (0.8) 0.4(0.1) 
0.15 (0.5) 0.2(0.1) 




Z{0)m 


Li(nH) 


L2(nH) 


CLci(fF) CiC2(fF) 


1-SQUID 
4-SQUID 


100 
100 


3.1 
3.8 


0.12 
0.3 


240 100 
240 24 



TABLE I: Parameter values used to fit the experimental 
I — V curves. The values resulting from independent measure- 
ments- are given in parentheses. We note that the experimen- 
tal value for Eji in the 1-SQUID experiment, obtained via a 
normal state resistance measurement, is unreliable because it 
was obtained after the probe was accidentally broken. When 
modelling the 1-SQUID sample, we used an asymmetry factor 
d — 0.12 for Eji, i. e. the Josephson coupling energies of the 
two JJs inside a SQUID satisfy \Ejji - Ejj2\/Eji = 0.12. 



Lorentzians to fluctuating lines. Therefore, we again as- 
sume that the measured I — V curves are a result of some 
kind of averaging over the quasicharge space, which now 
is two-dimensional because of two quasicharges Qq and 
Qq. The theoretical and experimental I — V curves are 
compared at <& = 0.45$o in the inset of Fig.|Sl One can 
see that they compare fairly well, even though the peak 
heights for show little disagreement. In this sample the 
band structure does not lead to strong peak splitting in 
contrast to the 1-SQUID sample. The reason for this is 
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FIG. 8: The current across the 4-SQUID sample when $ = 
(main frame) and $ = 0.45$o (inset). The solid line is the 
experimental and the dashed line the theoretical I — V curve. 
For the theoretical fit we have averaged over all values of 
Qo and Q'g and the widths of the bands are indicated as ar- 
rows. In the inset one can see that the |1, 0, 0)-resonance has 
widened and centers at ~ 32 nV. The |2*, 0, 0)-resonance lies 
almost at the same point as the upper edge of the |2,0, 0)- 
resonance, making the overall structure asymmetric. The 
|3*, 0, 0)-resonance lies between ~ 88 /iV and ~ 108 /iV, 
and the |3, 0, 0)-resonance lies between ~ 78 fiY and ~ 112 
fiV. In the calculation of the relaxation rates we have used 
Qint = 0.2(Qi + Qs + Qlci + Qlc2) (section|ryj. 



that the van Hove singularities in two dimensional qua- 
sicharge space are weaker than in one dimension. This 
is not the case for single excitations |1,0, 0), |2,0,0) and 
1 3, 0, 0), and indeed the band edges of the latter two are 
seen as separate peaks in the theoretical I—V curve. We 
have used uniform quasicharge distribution, which is the 
simplest guess as the physics of the average processing 
are unknown. 

The experimental and theoretical peak broadening are 
compared in Fig.|51for three of the transitions. The theo- 
retical width is obtained by summing up the peak width 
at <I> = 0, which for this choice of parameters is « 10 /^eV, 
and the increase due to broadening of the bands. There 
is good agreement between the theory and experiment. 
For example, the |2, 0, 0)-resonance broadens faster than 
the |2*, 0, 0)-resonance which is consistent with the theo- 
retical model. The width of the |3, 0, 0)-resonance is not 
analyzed since it rapidly becomes unobservable due to 
strong broadening. 



VI. CONCLUSION 

We have carried out a theoretical study of Cooper pair 
tunneling across a voltage biased asymmetric SCPT and 
a system consisting of three JJs in series, where the mid- 
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FIG. 9: Experimental linewidths (data points) compared with 
theoretical linewidths (lines) for three of the resonances as 
a function of the magnetic flux. The theoretical values are 
obtained by adding a constant term 10 /^eV to the widths of 
the bands. 



and band structure are taken into account. We pointed 
out, that the nature of the tunnehng across the probe 
turns from incoherent to coherent when the golden rule 
tunneling times exceed the relaxation times induced by 
the dissipative environment. Furthermore, we discussed 
that a simple master equation correction to the popula- 
tion of the eigenstates in the incoherent calculation leads 
to a good approximation for the current for arbitrary 
values of the voltage and for different flux values. 

In the last part of this paper we showed that a detailed 
theoretical understanding of experimental data can be 
achieved. In particular, the multiphoton processes be- 
tween different mesoscopic elements and spurious LC- 
resonators as well as the band structure of the Josephson 
junction can be probed by a small Josephson junction 
coupled to SQUID (s). Especially, the detection of en- 
ergy bands of higher excited states is confirmed by the 
fact that the observed widening of the resonances was in 
good accordance with the linewidths obtained from the 
model. 



die one acts as a probe, and applied the models in analyz- 
ing the experimental findings of Ref . [a. The treatment of 
the problem was done in the weak coupling regime, where 
the Cooper pairs tunnel incoherently across the probe, 
and was based on the idea of extending the well known 
P(£^)-theory into the regime where the anharmonicity 



Acknowledgments 

This work was financially supported by the National 
Graduate School in Material Physics, Finnish Academy 
of Science and Letters (Vilho, Yrjo and Kalle Vaisala 
Foundation) and the Academy of Finland (National Cen- 
ter of Excellence Programme). 



Electronic adress: |juha.leppakangas@oulu.5| 

G.-L. Ingold and Y. V. Nazarov, Single Charge Tunneling, 

ed. H. Grabert and M. H. Devoret (Plenum, New York, 

1992). 

G. L. Ingold, H. Grabert, and V. Eberhardt, Phys. Rev. B 

50, 395 (1994). 

E. B. Sonin, Phys. Rev. B 70, 140506 (2004). 

R. K. Lindell, J. Delahaye, M. A. Sillanpaa, T. T. Heikkila, 

E. B. Sonin, and P. Hakonen, Phys. Rev. Lett. 93, 197002 

(2004). 

R. Lindell, J. Penttila, M. Sillanpaa, and P. Hakonen, 

Phys. Rev. B 68, 052506 (2003). 

R. Lindell, J. Penttila, M. Paalanen, and P. Hakonen, 

Physica E 18, 13 (2003). 

R. Lindell Ph.d. thesis, Helsinki University of Technology, 

(2005). 

D. V. Averin, A. B. Zorin, and K. K. Likharev, Sov. Phys. 

JETP. 61, 407 (1985). 

K. K. Likharev and A. B. Zorin, J. Low Temp. Phys. 59, 



347 (1985). 

U. Weiss, Quantum Dissipative Systems, 2nd ed. (World 

Sci., Singapore, 1999). 

A. M. V. d. Brink, A. A. Odintsov, P. A. Bobbert, and 

G. Schon, Z. Phys. B 85, 459 (1991). 

A. M. V. d. Brink, G. Schon, and L. J. Geerligs, Phys. Rev. 

Lett. 67, 3030 (1991). 

T. Hoist, D. Esteve, C. Urbina, and M. H. Devoret, Phys. 

Rev. Lett. 73, 3455 (1994). 

J. Leppakangas and E. Thuneberg, to be published in AIP 

proceedings of LT-24. 

A. Shnirman and Yu. Makhlin, JETP. Lett. 78, 447 (2003). 

Yu. A. Pashkin, T. Yamamoto, O. Astafiev, Y. Nakamura, 

D. V. Averin, and J. S. Tsai, Nature 421, 823 (2003). 

O. Astafiev, Yu. A. Pashkin, Y. Nakamura, T. Yamamoto, 

and J. S. Tsai, Phys. Rev. Lett. 93, 267007 (2004). 

M. Tinkham, Introduction to Superconductivity, 2nd ed. 

(McGraw-Hill, New York, 1996). 



